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^ , Abstract 

The D3/D7 holographic model aims at a better approximation to QCD by adding 
£> ' to N = 4 SYM theory Nf of N = 2 supersymmetric hypermultiplets in the fun- 

damental representation of SU(N C ) - the "flavor fields" representing the quarks, 
ly-j ! Motivated by a recent observation of the importance of the Wess-Zumino-like (WZ) 

term for realizing the chiral magnetic effect within this model, we revisit the phase 
Q^ ■ diagram of the finite temperature, massless D3/D7 model in the presence of ex- 

ternal electric/magnetic fields and at finite chemical potential. We point out that 
the A-V-V triangle anomaly represented by the WZ term in the D7 brane probe 
action implies the existence of new phases that have been overlooked in the previous 
studies. In the case of an external magnetic field and at finite chemical potential, 
rN \ we find a "chiral helix" phase in which the U(1)a angle of D7 brane embedding 

increases monotonically along the direction of the magnetic field - this is a geo- 
metric realization of the chiral spiral phase in QCD. We also show that in the case 
of parallel electric and magnetic fields (E, B) there exists a phase in which the D7 
brane spontaneously begins to rotate, so that the U(1)a angle changes as a function 
of time - this may be called the "spontaneous rotation" phase; it is a geometrical 
realization of a phase with non-zero chiral chemical potential. Our results call for 
a more thorough study of the (T, B, E, fi) phase diagram of the massless D3/D7 
model taking a complete account of the WZ term. We also speculate about the 
possible phase diagram in the massive case. 
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1 Introduction 



This work is motivated in part by a recent realization of the chiral magnetic effect (CME) 
[U [2] in the D3/D7 model via the Wess-Zumino-like (WZ) term in the D7 probe brane 
action [3]*]. The field theory dual to the D3/D7 holographic model is the N=4 super- 
symmetric Yang-Mills (SYM) theory with Nf of N=2 hypermultiplets in the fundamental 
representation of the color group SU(N C ) - the "flavor fields" representing the quarks. In 
this dual field theory, the axial phase of a Dirac quark in the hypermultiplet is tied to a 
U(1)r subgroup of the SO(6)r symmetry of the N=4 SYM sector. The hypermultiplet 
also has its own flavor U(l)v symmetry which is vector-like to the hypermultiplet quarks. 
These ingredients are sufficient to give rise to the U(l)yU(l)n flavor anomaly through 
the triangle one-loop diagram involving the hypermultiplet quarks. In the holographic 
picture of D7 probe brane in the large N c D3 brane background, this triangle anomaly is 
captured by a WZ-like coupling in the D7 probe brane action. A peculiar feature in this 
set-up is that £7(1) r (we will also call it U(1)a interchangeably) is geometrically realized 
as a £7(1) angle in the 5- sphere of AdS^ x S 5 geometry, whereas the dynamics of £7(1) y 
is as usual described by the gauge field on the D7 brane - therefore the triangle anomaly 
is no longer represented by a 5D Chern-Simons term. 

Because £7(1)# is shared by other adjoint matter fields in N=4 SYM, the U(1)r charges 
in the hypermultiplet can be lost to the adjoint sector, and the axial chemical potential 
can be meaningful only in a quasi-equilibrium sense. A novel idea in Ref. j3] is to introduce 
an external time-dependent axial phase into the system that simulates an axial chemical 
potential. As U(1)a is geometrically realized as an angle in the 5-sphere, this can be 
achieved by externally rotating the D7 brane along the £7(1) a angle. Considering that 
the axial charge can be lost to the adjoint sector, it is a nice way of maintaining a chiral 
imbalance in an otherwise equilibrated system [3] . This is in fact similar to the true axial 
£7(1)a symmetry of real QCD where at finite temperature the sphalerons can create or 
annihilate the axial charge. Indeed, a time- dependent QCD #-angle which is equivalent to 
a time-dependent axial phase by anomaly was one way of effectively describing the axial 
chemical potential, with /jla — j9]. Since in the vacuum picture 6 can be interpreted 
as a quasi-momentum (in analogy to the Bloch crystal), 9 may be interpreted as an 



* Closely related phenomena have been discussed in the physics of neutrino emission [3], primordial 
electroweak plasma [5] and quantum wires [6]; the separation of electric charge in QCD plasma induced 
by the chirality imbalance in the presence of magnetic field and/or angular momentum was first discussed 
in 018]. 
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external "force" that maintains an imbalance between the left and right fermions. This 
can also be realized holographically, for example in the Sakai-Sugimoto model, through a 
time-dependent Cjyj holographic dual to Qqcd which could clarify the treatment of the 
holographic chiral magnetic effect [10J QH OH H31 El H51 M\ ■ 

Indeed, one issue with the holographic description of the CME is that in Minkowski 
signature black hole space-time, the temporal component of the vector potential at 
the UV boundary is not quite equivalent to the chemical potential even in the case of 
non-anomalous global symmetries [T4"| ITT] . The thermal Green's functions in Minkowski 
signature black hole space-time have the structure 

Tr ensemble(C , (t)C , (i / )) time ordering > C 1 - 1 ) 

where 0(t) is an operator that evolves as 

Q(t) = e +iAt O(0)e- iAt , (1.2) 

while the ensemble trace can be either grand canonical or (micro) canonical, 

Tr enscmble = Tr (e-^-^) or Tr (e"^) , (1.3) 

where N is the charge operator. A crucial point is that the hamiltonian appears here in two 
places with rather different roles - the one appearing in the ensemble trace, e~^ H ~^ N ^ or 
e~^ H , is responsible for the Euclidean imaginary time evolution of period j3 while the one 
appearing in 0(t) is responsible for the evolution in Minkowski real time. In Minkowski 
black hole space-time, one sees only the Minkowski evolution and the Euclidean one is 
not manifest geometrically. Conversely, in the Euclidean black hole geometry, only the 
Euclidean evolution is geometrically realized. To be more explicit, by turning on A t (oo) at 
the UV boundary in Minkowski signature black hole, one achieves a famous replacement 
according to the AdS/CFT dictionary, 

H^H-A t (oo)N ; (1.4) 

the question however is which hamiltonian should we replace on the field theory side: the 
one in the ensemble average or the one in 0(t), or both? Since the boundary field theory 
is in Minkowski signature space-time, it seems clear that the hamiltonian in 0(t) has to 
be replaced; however the question whether the hamiltonian in the ensemble average also 
has to be replaced is more difficult to answer. To establish the identification A t (oo) = fi 
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in the ensemble average one should go to the Euclidean signature black hole where the 
imaginary time evolution becomes geometrically realized. 

It is clear from the above discussion that one should not naively associate A t (oo) in 
Minkowski black hole space-time with a chemical potential. This is relevant because the 
ambiguities of the chiral magnetic current appear precisely due to A t (oo) in Minkowski 
signature space-time, and are likely related to its effect on 0(t) which should not have 
been the case for the true chemical potential. Without A t (oo) the ambiguities of the 
holographic CME disappear [HI [17]. In Minkowski black hole, the value of the chemical 
potential can be obtained by the work done to a unit charge in bringing it from the UV 
boundary to the black hole horizon. The issue of singularity at the horizon poses no 
problem because what matters in Minkowski dynamics is only the future event horizon 
and the gauge field is regular there for any choice of A t (oo) [17J. This is in accord with 
the fact that A t (oo) in Minkowski black hole is not quite the chemical potential. One 
may choose to perform a bulk gauge transformation to remove A t while introducing a 
time-dependent A r instead, 

A r = -td r A t (r) . (1.5) 

Physically, this introduces a time-dependent Wilson line stretching from the UV boundary 
to the horizon, 

W = e i ^B Ar = e tt ^ , (1.6) 

with frequency which is equal to the chemical potential. What is important here is that 
there is no singularity issue with this Wilson line. In fact, the norm AmAnq mn one uses 
in the singularity argument is not a gauge- invariant concept. 

Although one needs an axial chemical potential to observe the CME, there are other 
interesting phenomena stemming from the triangle anomaly. One example is the chiral 
separation effect [TBI EE]: the emergence of an axial current along the magnetic field 
B = Bx 3 in the presence of an "ordinary" vector chemical potential ^y, 

A = §t»vB • (1.7) 

The close connection between the CME and the chiral separation effect is particularly easy 
to see in the case of a strong magnetic field when dimensional reduction is appropriate 
[2D] . In the dimensionally reduced (1 + 1) theory, the axial and vector currents are related 
by jfi — e ^vj% so that the axial charge density j Q A induces a vector current j 1 (CME) and 
the vector charge density j° induces an axial current j\ (chiral separation). 
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Contrary to the axial chemical potential, the vector chemical potential (as well as 
magnetic field) is much easier to introduce in the D7 brane. Once we accept (11.7j) . it 
is natural to look for a signal of the axial current j\ in the following way: as £7(1) a is 
geometrically realized the axial current will take a form of a chiral spiral [2TJ [22j [23j EU |25j 
EH EE] that is, the £7(1) a angle of the D7 brane embedding shape should have a constant 
gradient along the space direction x 3 = z. Even though one is working in the massless 
limit, the D7 brane should develop a profile of non-zero axial phase gradient along x 3 to 
satisfy the constraint (11. 7p dictated by the triangle anomaly. We will call this the chiral 
helix phase emphasizing its geometrical realization in the holographic setup. 

The phase diagram of D7 brane dynamics in the presence of both magnetic field and 
ordinary chemical potential at finite temperature has been studied before (27J EH] , but 
the triangle anomaly constraint (II. 7p . or equivalently the effect of the WZ term in the 
D7 brane action, seems to have been missed in these analyses. In this paper, we find a 
significant modification of the phase diagram due to the anomaly, and prove the emergence 
of the chiral helix phase through a dynamical instability. A more thorough study of the 
full {B,n) phase diagram will be presented elsewhere [29] . 

Another interesting and well-known effect stemming from the triangle anomaly is 
the creation or annihilation of the £7(1)^ charge in the presence of parallel electric and 
magnetic fields, 



At finite chiral chemical potential /i^, ( 11.81) determines the power of the chiral magnetic 
current in the case of the parallel electric and magnetic fields: 



Note that no power is dissipated in the absence of the electric field, and that depending 
on the relative signs of E ■ B and /j,a, the power can be either positive or negative - in the 
latter case the current is powered by the energy stored in the system due to the difference 
of Fermi energies of left and right chiral charges [21 [9] . The reversibility of the sign of the 
power signals the lack of dissipation that is a consequence of the time reversal invariance 
of chiral magnetic conductivity and other anomalous transport coefficients - this provides 
an important constraint on the anomalous hydrodynamics [30] . 

In the holographic setup one can easily introduce the electric/magnetic fields through 
world- volume gauge field on the D7 brane. Because the £7(1) a charge is simply the 
angular momentum of the D7 brane along the axial angle in the 5-sphere, the above 



d^-E-B^O 



(1.8) 




(1.9) 
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triangle anomaly constraint implies that the D7 brane should start rotating along this 
U(1)a angle with an increasing speed. Due to the loss of axial charge (dissipation) to 
the adjoint sector, one expects the system to stabilize with a finite angular momentum 
afterwards. We will call this the "spontaneous rotation phase". The existence of this 
phase seems to have been missed in the previous analyses in the literature [27J IS]; we 
expect that the spontaneous rotation phase leads to a significant modification of the phase 
diagram. In the present work we prove the existence of the spontaneous rotation phase 
by observing a dynamical instability towards it. The study of the full phase diagram of 
(B, E, fi) is deferred to future. 

We stress that there are no externally driven time-dependent parameters in our situ- 
ations, contrary to Ref. [3]. This is because the chemical potential and electric/magnetic 
fields of the baryonic U(l)y symmetry are easier and more natural to introduce through 
the D7 world- volume U (l)y gauge field. Yet we observe that the existence of the U (l)yU (1)r 
triangle anomaly represented by the WZ term still leads to interesting consequences. 

Although we focus only on the massless case in this paper, it is interesting to speculate 
about the massive case. In the case of the chiral spiral in the chiral-symmetry broken 
phase of QCD, the massive case features an inhomogeneous chiral spiral : the axial phase 
jumps only in narrow periodic ranges of x 3 [22]. This is due to a competition of the free 
energy associated with the phase gradient and the anomaly constraint imposed by (|1.7j) . 
One naturally expects that a similar phenomenon would happen in the massive D3/D7 
model, and it will be interesting to pursue this further. 



2 Massless D3/D7 model with (T, B, fi) 

We study the dynamics of the probe D7 brane embedded in a gravity background of large 
N c D3 branes at finite temperature T given by 



ds 2 = ^ (-V(r)df- + YXdx'yJ + ~^rf 2 + £ 2 <«i I , V(r) = 1 - ^ 

where L A = 4-Kg s N c lj = A/ 4 and €5 is the volume form of a unit 5-sphere; g s and l s are 
the string coupling and string length. In the weak coupling limit, the D3 branes span 
(t,x l ),i = 1,2,3 and the D7 brane wraps additional four dimensions we call x 4,5 ' 6,7 . The 
remaining two dimensions x 8,9 are transverse to both D3 and D7 branes. The six dimen- 



5 




Figure 1: A schematic picture of D7 brane embedding in the D3 background. The (x 8 , x 9 ) 
plane of our interest is shown explicitly. 

sional space x A ~ 9 is the total transverse space to the D3 brane, and the radial direction r 
and the 5-sphere O5 in (I2.10p are the radius and angles in this space. The theory of the 
D3/D7 branes is N=4 SYM theory of SU(N C ) plus N=2 hypermultiplet of fundamental 
representation fields. In the strong coupling regime that we are interested in, the D3 
branes are replaced by the above gravity- flux background ( 12.101) whereas the D7 brane is 
treated as a probe brane to this holographic background. On the field theory side, the 
dynamics of the probe D7 branes should be dual to the dynamics of N=2 hypermultiplet 
at strong coupling in the quenched approximation. 

The embedding geometry of the D7 probe brane in the above background ( 12. 10}) is 
explained in FigJTJ In general, the D7 brane bends in the transverse direction x 8,9 , and 
this is all the data one has to specify for the shape assuming that the shape is invariant 
under rotations in x 4 ' 5 ' 6 ' 7 space. Introducing the radius of x 4,5 ' 6,7 space that the D7 brane 
spans 



the D7 brane wraps a 3-sphere of constant p for each p. Note that p = p(r, x^) is in general 
a non-trivial function of holographic 5-dimensions (r, x^) depending on the bending in x 8,9 




(2.11) 
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directions by the relation 

p = yV 2 - (A 8 (r, x»)) 2 - (X 9 (r, x^)) 2 , (2.12) 

where X 8,9 (r, x^) are functions that specify the bending shape in x 8 ' 9 space. One can con- 
veniently choose the eight-dimensional D7 brane world- volume coordinates to be (x M , r, f^), 
and the embedding geometry is completely determined by two functions X 8,9 (r, x^). 

The rotational U(l) in x 8,9 space is a part of 5'0(6)-symmetry of the total transverse 
space x 4-9 which is the holographic manifestation of SO(6)r R-symmetry of N=4 SYM 
theory. For the hypermultiplets it also corresponds to a phase rotation of the N=2 super- 
symmetric complex mass term, as the r — > oo asymptotic value of (X 8 + iX 9 ) is precisely 
such mass parameter for the hyp er mult iplet. Although the hypermultiplet involves scalars 
too, at least for the (Dirac) fermions it is very similar to the axial U(1)a in real QCD. We 
therefore call it either U(1)r or U(1)a interchangeably. The novelty here is that U(1)a 
is geometrically realized as a real rotation in x 8 ' 9 internal space. Several previous stud- 
ies have explored this aspect to get useful results relevant for chiral symmetry breaking 
phase transitions in QCD [221 E2 EH ES] B More relevant to our present work, Ref.[3] 
recently simulated the axial chemical potential in QCD by introducing an external time- 
dependent rotation along this x 8 ' 9 U(l) angle. One drawback is that this R-symmetry 
is shared by adjoint scalars and fermions in N=4 SYM theory, so that the total charge 
can be lost to the background geometry |3j. However, the effective chemical potential 
that they introduce by rotation can be kept stationary by continuous external inflow of 
necessary charges and one can meaningfully discuss the physics of hypermultiplet sector 
with a finite R-symmetry chemical potential. 

This bears some similarity to real QCD where axial charges may be lost due to QCD 
sphalerons, and indeed a time-dependent external 6*q<7£t angle (which is equivalent to the 
axial phase by anomaly) was one of the ways to introduce an effective axial chemical 
potential [9]. It is perhaps relevant to point out that one can do a similar thing even 
in the Sakai-Sugimoto model. The U(1)a suffers non-conservation due to coupling to 
the C RR which is dual to 0Qc#-angle. This is a holographic manifestation of U(1)a 
anomaly with QCD gluons. Although this effect has been neglected based on large N c - 
suppression, one can go on to introduce a time-dependent C RR to introduce an effective 
axial chemical potential; this procedure would bypass the issues regarding the holographic 
chiral magnetic effect discussed in the Introduction. 

tAs in QCD, U(1)a is not a true symmetry due to the anomaly involving gluons. A typical justification 
for discussing chiral symmetry breaking with it is a large N c suppression of this anomaly. 
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The D7 brane dynamics also includes a baryonic U(l) symmetry by a U(l ) gauge field 
residing on its world-volume. The situation we are going to study is the one having a 
constant magnetic field along, say x 3 = z direction, 

F 12 = B , (2.13) 

and a finite chemical potential /U, both of which are with respect to this baryonic symmetry. 
Note that B which is constant everywhere is a trivial solution of equations of motion. We 
will focus on the massless case which means 

(X 8 ,X 9 )^0 , r->oo , (2.14) 

and we will discuss what might be happening in the massive case in the last section. The 
phase diagram of the system in (T,B,fj) has been studied previously in Ref . [28] . but as 
discussed in the Introduction an interesting role played by U{1) 2 U(1)r triangle anomaly 
represented by a Wess-Zumino term was overlooked. Given (B,[x), the triangle anomaly 
dictates an existence of the chiral separation effect: 

f R = ^B , (2.15) 

and in our case of U(1)r geometrically realized as a rotation in x 8 ' 9 space, a finite j\ 
current would take a form of a helix, i.e. a non-zero spatial ^-gradient of the U(1)r angle 
of the D7 brane embedding, see Figj2] for an illustration. This is an analog of the chiral 
spiral of pion gradient in low-energy QCD with electromagnetic B and \Ib [22], except 
that the axial phase and the spiral in our case are realized geometrically and are easily 
visualized. 

This implies that a sizeable fraction of the phase diagram with (B, //) should in fact be 
the chiral helix phase. To establish its location in the (T, B, /i) phase diagram, one should 
compare the grand canonical free energies (including the effects from the WZ term) of the 
phases with and without the chiral helix. This interesting task will be pursued elsewhere 
[29] . but in this paper we will prove the existence of this phase by showing a dynamical 
instability towards it from the phase which does not possess the chiral helix. FigJH] is a 
rough picture of the phase diagram of (T,B,/i) in Ref. [28] without considering the WZ 
term. Due to conformal symmetry, the only meaningful parameters are (-^g, ^) arid one 
The vertical axis is related to T by 

T2^ T ■ (2M) 



sets B = 1 



* Precisely speaking there is also a rescaling (2itl 2 s )F — > F and we choose l s to have L 4 = AZf 
before setting B = 1. We also will use this convention later. 
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Figure 2: The D7 brane embedding in the chiral helix phase. The U (1)r angle 4> in (x 8 , x 9 ) 
plane is monotonically increasing along x 3 = z direction. In the space of (x 8 ,x 9 ,z) the 
shape indeed looks like a helix. 

The region III in the upper-right part is the so-called supersymmetric embedding phase 
where the D7 brane embedding is straight in x 8 ' 9 space without bending all the way to 
the black hole horizon, (X 8 (r), X 9 (r)) = (0,0). We are going to study the dynamical 
instability of this phase towards forming the chiral helix, and as an exemplar case we will 
pick one point - the point A in Figj3]- which has (^|,^) = (0.25, 1) and is well inside 
the region III without any ambiguity. We will indeed find that this point is unstable to 
linearized chiral helix modes. 

The action of D7 brane is 

S D7 = -fPr J A e-ydet(& + 2 nP s F) + Vi^f- j Cf R A F A F , (2.17) 

where /x 7 = (27r)~ 7 Z,7 8 , = g s , and F is the field strength of the baryonic U(l) gauge 
field. Note that the embedding dynamics of X 8,9 (r, x M ) enters through the induced world- 
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Figure 3: A schematic picture of (T, B, fj,) phase diagram without considering triangle 
anomaly (WZ term) from Ref.[28| [31]. B has been set to 1. The point A has (^ ; ^) = 
(0.25, 1) and has an instability toward chiral helix phase. 

volume metric g*. Because l s drops in any field theory observables, one can conveniently 
choose it such that L A = A/ 4 = 1, and we also rescale {2itl 2 s )F — > F for simplicity. We are 
going to study linearized fluctuations from the configuration of point A in Figj3] which 
has X 8 ' 9 (r, x^) = 0. The chemical potential, or equivalently the background solution of 
F tr is obtained from the action 

S D 7 = -^f / d 5 x {rV^+sV 1 - (F tr ) 2 } , (2.18) 
which gives one in A r = gauge, 

F<?> = -M 0) = —, Q , (2-19) 

lQ 2 + r 2 (r A + B 2 ) 



where a constant of motion Q is determined from fi by the condition 

H = Af ] (oo) = r dr Q = . (2.20) 



y/Q 2 + r 2 (r 4 + B 2 ) 

It is tedious but straightforward to expand the action (I2.17P quadratically from the back- 
ground solution in terms of small linearized perturbations of the modes (X 8 ' 9 , A t , A z ) (we 
omit 5-symbol for simplicity) assuming the dependence only on (t, z, r) that are poten- 
tially relevant for the chiral helix instability. We have verified the consistency of this 
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ansatz. We find that gauge field perturbations decouple from those of X 8 ' 9 so from now 
on we keep only X 8,9 perturbations that are of interest for us. One subtlety regarding 
C RR should be mentioned: writing the 5-sphere metric as 



dflj 



dx 2 + sin 2 xd<fi 2 + cos 2 xd^l 



where < x < § is the angle from the p axis such that (see Fig{T]) 



P 

cos x = - 
r 



R 

sin x = — , 

r 



R 



^(x 8 ) 2 + (x 9 ) 2 



and (p is the U(1)r angle in x 8,9 -plane, the volume form takes a form 

1 



£5 = sin x cos x dx A d(p A £3 = d 



4 



cos x A d(j) A e 3 



where 63 is the volume form of unit 3-sphere. From this one obtains C RR as 



[2ixl s ) A N c 1 



(C - cos 4 x) A d<p A e 3 , 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



7p3 4 

where a constant C is a freedom of sigular gauge transformations. Since the background 
D7 brane shape X 8,9 (r) = corresponds to x — and we are looking at small fluctuations 
around it, we need to choose C such that C RR is regular around x — 0- At x = 0, the 
angle </> becomes singular, and Cj^ should vanish to be regular, which fixes C — 1. (The 
other choice C = would make ^ regular at % = | instead where the 3-sphere vanishes. 
This will be suitable when discussing fluctuations around Minkowski embeddings.) This 
finally gives us an expression for the WZ term: 
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wz 



1287T 4 



d 5 x 



2r 2 - (X 8 ) 2 - (X 



MNPQR ( X S 9mX 9 _ X 9 d M X 8 ) F NP F QR , 

(2.25) 



where ^ MNp Q R is purely numerical. 

After a sizable amount of algebra, the quadratic action for the fluctuations one gets is 
XN r r , K fl.,,/„„ B \2 1„_/„„ n2 „ x2 „ /(A 8 ) 2 \ 



167T' 



| d 5 x \ l -A{r) A 8 ) 2 - is(r) (<9 2 X 8 ) 2 - l -C{r) ((d r X 8 ) 2 - d r 



+ ^D(r)(A 8 ) 2 + (same with X 8 ^ X 9 ) - ^E(r) [X 8 8 Z X 9 - X 9 8 Z X 8 ) j , (2.26) 



where the coefficient functions are given by 

1 1 



A(r) 



Vr 4 + B 2 



r A V(r) 



F 



(0) 



tr 



. B(r) = rVr A + B 2 — Jl 



7(0) 

tr 
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tr 



C{r) = rVr* + B 2 l \ J , D(r) = r^r 4 + B 2 — \jl - [F t 

E(r) = ^4 0) . (2.27) 

The last term results from the WZ term which plays an essential role in our discussion. 

It is a straightforward exercise to derive the equations of motion from the above and 
to study the linearized stability of the system. As one has a black hole horizon located at 
r = th, one needs to impose a physically relevant in-coming boundary condition at the 
horizon, and for this purpose it is more convenient to work in the Eddington-Finkelstein 
coordinate (t*,r*), 

f r dr' 

U=t + 7 /\ o t t / ,\ , n=r , (2.28) 



oo (r') 2 V(r') 



upon which one has 



d t = d u , d r = d rt + „} rl d u , / dtdr = I dt^dr* . (2.29) 



r 2 V (t 

The quadratic action (I2.26P then takes a form in Eddington-Finkelstein coordinate (we 
omit subscript * for simplicity), 

- l - Q;(d r C(r)) - D(r)) (X 8 ) 2 + (same with X 8 <-> X 9 ) 

- l -E{r) (x*d z X 9 - X 9 d z X s ) | , (2.30) 
Note that (d t X 8 ' 9 ) 2 terms disappear completely due to the identity 

A ^-7WW • (2 ' 31) 

and all the coefficient functions that appear in the above, especially J^y^ , are regular at 
the horizon r = r#. 

Recall that the usefulness of the Eddington-Finkelstein coordinate is that simple regu- 
larity at the horizon r = m automatically guarantees the in-coming boundary condition, 
while out-going modes look singular in the coordinate. An easier way to see this is to 
consider a wave oscillating in time t* in the Eddington-Finkelstein coordinate 



) , (2.32) 
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which automatically contains the necessary in-coming radial phase part in terms of the 
original Schwarzschild coordinate, so that any regular wave in Eddington-Finkelstein co- 
ordinate is in-coming at the horizon. 

Equations of motion from (I2.30P are 

B(r)#t*» + d r {C{r)d r X*») + ^OAX* + d t d r (^^") 
- (- (d r C(r)) - D(r)) X 8 ' 9 t E(r)8 z X 9 ' 8 = , (2.33) 



where the last term mixes X 8 and X 9 in such a way that the correct diagonal basis is in 
fact a helical basis = X 8 ± iX 9 , which is also tied to the z-direction flipping as one 
expects. One can easily see that the modes at hand indeed correspond to the chiral helix 
modes with finite ^-momentum. One needs to consider only X^ + ' = X in the analysis as 
the equations are invariant under X^ — > X^ and z — > —z. 

To study instability associated with a finite ^-momentum one goes to frequency- 
momentum space by assuming e ~ luJt + lkz ^ after which one has 

C{r)dlX + [{d r C{r)) ~ 2 ^^)) 

- (^-(d r C{r))- D{r) + k 2 B(r) + kE(r) + iud r ^-^LjW=0 . (2.34) 

As mentioned before, one imposes regularity at the horizon r = r# for in-coming boundary 
condition. Note that C{th) — and the boundary condition is not trivial. On the UV 
boundary r — > oo, one can have two possible modes from the above equation as usual in 
AdS/CFT, 

X ~ X - iu^- + ^ + • ■ ■ , (2.35) 

where X corresponds to the bare mass of the hypermultiplet and the normalizability con- 
dition puts the constraint X = in our massless case. These two boundary conditions 
give a discrete spectrum of u for a given k, which is an example of quasi-normal mode 
problem. If the lowest u(k) — > in the limit k — > 0, this is sometimes called hydrody- 
namic dispersion relation for an obvious reason, but here we are interested in the finite k 
spectrum. An inspection of our master equation f)2.34p shows that u is purely imaginary, 
so that one can write 

u) = ilm[u} . (2.36) 
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Considering e~ tuJt dependence, any positive Im[u;] > signals an exponentially growing 
amplitude and hence an instability, while one expects Im[oj] < for typical dissipative 
relaxations. 

Before delving into numerical searches for instability that will be described shortly, 
one can qualitatively understand from (I2.34p how the unstable modes of Im[u;] > can 
possibly appear. Consider first a fictitious situation where there was no WZ term, or 
equivalently let E(r) be absent in (12.341) . This would be a typical situation of D7 brane 
world-volume fluctuations along transverse X 8 ' 9 directions, and one can expect that the 
modes would simply be dissipating in the presence of black hole. The dependence on k 
should obviously be that the larger k 2 is the faster the mode decays, so that the quasi- 
normal spectrum would be 

lm[u] = -miff -ak 2 H , a > , (2.37) 

for a reasonable range of k. Then, in (12.341) the effect of having WZ contribution is simply 
replacing B(r e ff)k 2 with B(r e ff)k 2 + E(r e ff)k where r e ff is the most relevant value of r 
for the mode wave function. Therefore the quasi-normal mode spectrum including WZ 
contribution should qualitatively look like 

lm[u] = -m 2 eff - ak 2 - (3k + ■ ■ ■ , a>0 , (2.38) 

where (3 is roughly proportional to E(r e ff). Completing square of the expression, one 
arrives at 



lm[u] = -a(k + ^-\ +^-m 2 eff = -a(k-k ) 2 + ^-m 2 eff , (2.39) 



so that when ^ — m\ff > 0, Im[u;] can be positive for a range of non-zero k centered 
around k = k$. This qualitative picture is confirmed in our numerical studies. 

In our numerical simulation, we solve (I2.34p from r = r# + £ to a UV cutoff r = r max 
given (&,Im[a>]) using a Mathematica package NDSolve. Note that regularity at r = r# 
from <^M) fixes the ratios and at r = th unambiguously, which allows one 

to start from r = th + s with a small number e. As the equation is linear we are free 
to choose the normalization X(th) = 1. We then impose a condition of normalizability 
|X(r ma:r )| < e' with another small number s' . We choose the parameters 

£ = 0.001 , r max = 10 , e = 0.05 , (2.40) 

in our numerical result in FigJH which clearly shows a range of k with positive Im[u;]. 
This is a numerical proof of the existence of the chiral helix phase. 
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lm(a>) 




-0.4 - .* 



-0.6 u 

Figure 4: Our numerical result of Im[u;] as a function of z-momentum k. We set B = 1 
and (t^jA 4 ) = (0.25,1). This shows that Im[w] > for a range of k and indicates an 
instability. 

3 Massless D3/D7 model with (T, £, £) 

In this section, we study another case of interest that should also be affected by the 
triangle anomaly, or equivalently the WZ term: the case when both electric and magnetic 
fields are present, and parallel to each other. The triangle anomaly dictates that the axial 
U(1)r current is not conserved under this situation 

d,f R ~E-B , (3.41) 

and one expects a continuous creation (or annihilation) of U(1)r charges in the system. 
Noting that the U(1)r charge at hand should correspond to the angular momentum of the 
D7 brane along the x 8 ' 9 -plane, the only way to satisfy this anomaly constraint is to have a 
spontaneous rotation of D7 brane in x 8,9 plane. Initially, the D7 brane angular momentum 
would increase; then, due to the existence of black hole horizon, the dissipation of the 
created U(1)r angular momentum to the background geometry (or adjoint sector) turns 
on, and one can imagine a stationary situation of a constant angular momentum carried 
by the rotating D7 brane. We will call this a "spontaneous rotation phase". Previous 
studies on the system [27J EE! seem to have missed this possibility, and a more complete 
study is certainly desirable [29J. In this section, we will prove the existence of this phase 
by showing a dynamical instability toward it from the configuration that does not possess 
rotation. 

The phase diagram of (T, B, E) from Ref.[3T] is sketched in Figj5] We focus only on 
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Figure 5: A rough picture of phase diagram of (T, B, E) without considering triangle 
anomaly in Ref. |31j. B is set to 1. The point A has (^,E) = (0.25,1), which will be 
shown to be unstable toward spontaneous rotation phase. 

the zero chemical potential case for simplicity, but the essential feature of the analysis 
in this section is independent of the presence of chemical potential. Using conformal 
symmetry, B has been set to unity, and the conventions are as explained in the previous 
section. The region III in the far right is again the supersymmetric embedding phase 
where the D7 brane is straight with X 8,9 (r) = 0. We will study the spontaneous axial 
rotation instability of the point A with (-^|, E) = (0.25, 1) which is well inside this region. 

To find the zeroth order background solution with constant (B, E) in supersymmetric 
embedding phase, one starts with the relevant action component 



Sm = ~T^J d " X {^^5^1 + V{r){F zr f - ^£ 2 } , (3-42) 

where one turns on external magnetic F12 = B and electric F tz = E fields along z = x 3 
direction. One can check that the constant B and E solve full equations of motion 
consistently, and one only needs to solve for F zr , or A z in A r = gauge. Solving it from 
the action, one obtains 



jr(o) J 



1 r 4 V(r)^ 2 



V(r)^(r* + ^)-^j 



(3.43) 



where J is an integration constant which is directly proportional to the induced baryonic 
current along z direction. One notices that the expression inside the square root in (I3.43P 
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can change signs if J is not suitably chosen [36]. The on-shell action itself also involves 
similar factors 



'on— shell 



XNc 

167T 4 



d b x { r 2 (r 4 + B 



1 ~ TW^PjE 2 
^ r2(r 4 + 52) __JLj ' 



(3.44) 



and one should make sure that the total expression inside square root remains positive 
for a meaningful solution. The numerator changes sign at r = r* where 



rt = r 4 H + E 



(3.45) 



and this should also be the point where the denominator changes sign, which determines 
J as 1 36 1 



E 1 



r 4 H + E 2 + 



B 1 



H 



+ E 2 



(3.46) 



The point r = r* will play an important role when we discuss linearized fluctuations from 
this background solution. 

Our task is to expand the D7 brane action quadratically for linearized fluctuations 
from this background solution. After some algebra, one finally arrives at 

S {2) = ^Jd 5 x t~A(r) (d t X 8 f + B(r) (d t X 8 ) (d r X 8 ) - l -C{r) (d r X> 

- i 0- (d r C(r)) - D(rj\ (X 8 ) 2 + (same with X 8 <-> X 9 ) 

+ ±E(r) (x 8 d r X 9 - X 9 d r X 8 ) + l -F(r) (x 8 d t X 9 - X 9 d t X 8 ) j , (3.47) 

where the last line is from the WZ term. The coefficient functions are 



A(r 
B(r 
C{r 

D(r 
E(r 



Vr 4 + B 2 



1 f 1 



r 4 \V(r) 



1 



r 4 V(r)' 



■E 



rv / r 4 + B 2 —F^ 



l + V(r)(F^) ~^E 2 



Vr 4 + B 2 \ (r*V(r) - E' 



1 + V(r) [F: 



\r-,E 2 



r 4 V(r) 



Wr 4 + B 2 - X 1 + V(r) [F s 



' zr 



r A V(r) 



E 2 



AB 



E 



F{r ) = l| F (o) 



„2 zr 



(3.48) 
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Figure 6: The shape of the D7 brane in the spontaneous rotation phase. The shape is 
helical along the radial direction while the system rotates in time. 



One should note that the expression 



k + V{r)(F^ ] ) 2 - - E 2 = ry/ r 4 + B 2 

^ > r*V(r) V \\ r 2 (r* + B 2 ) - ^ 



1 

r 4 V(r) 



E 2 



(3.49) 



is regular and non-vanishing at r = r* precisely because of the choice of J above, so that 
all coefficient functions are regular at r = r*. 

It is an important fact however that C{r*) = 0, which implies that the regularity 
boundary condition at r = r* is a non-trivial one. Combined with the UV normalizabil- 
ity boundary condition, these two boundary conditions give us a discrete quasi-normal 
spectrum of the fluctuations. Therefore r = r* (recall r* > r#), instead of the black 
hole horizon, seems to play a role of the IR boundary for linearized fluctuations on the 
D7 brane [37]. Because A(r) remains finite at r = r*, one does not need to work in the 
Eddington-Finkelstein-like coordinate around r = r*, and a regular wave at r = r* is not 
necessarily in-coming. 

The equations of motion are 
-A(r)d 2 t X^ - B(r)d t d r X 8 > 9 - d r (B(r)d t X^) + d r (C(r)<9 r X 8 ' 9 ) (3.50) 
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Figure 7: Numerical result for the lowest complex frequency of a linearized spontaneous 
rotation mode. The background has (^,E,B) = (0.25,1,1) in the supersymmetric 
embedding phase. This is a numerical proof of the existence of spontaneous rotation 
phase. 

i (d r C(r)) - D(r)j X 8 ' 9 ± ^E(r)d r X 9 > 8 ± l -d r (#(r)X 9 ' 8 ) ± F(r)d t X 9 > 8 = , 

and in terms of helicity basis X^ = X 8 ± iX 9 , the equations of motion become diagonal. 
This indicates that the D7 brane shape will be helical in the radial r direction, while 
rotating in time, see Figj6j To study the instability, we go to the frequency space assuming 



- iuit 



and find that under X^ <H- X^ ) the frequency maps to w H —to so that the 
imaginary part of u is the same for both JW. Since the instability is signaled by a 
positive imaginary part Im[o>] > 0, one needs to study X^ = X mode only. The 
equation to solve in frequency space is then 



C{r)d 2 r X + (<9 r C(r)) + 2iuB(r) - iE(r) )d r X 



(3.51) 



+ ( u 2 A(r) + icu (d r B(r)) - - (d r C(r)) + D(r) - - (d r E(r)) - uF(r) ) X = . 

r 2 



We solve (13.511) numerically from r = r^ + e until a UV cutoff r max given a complex uj. 
One can normalize X(r*) = 1 and the regularity at r = r* fixes the initial conditions. We 
then test a UV normalizability condition |V(r max )| < e', and if the solution satisfies it we 
let the program put a dot in the complex w-plane. FigJT] is our result using parameters 

e = 0.001 , r max = 10 , e' = 0.05 , (3.52) 



19 



which clearly proves that the quasi normal spectrum u has a positive imaginary part, and 
hence instability. 

4 Summary 

To summarize, we have studied the phase diagram of the D3/D7 holographic model in two 
cases: 1) at finite (vector) chemical potential and in the presence of an external magnetic 
field, and 2) in the presence of external parallel electric and magnetic fields. We have 
found that the triangle anomaly represented by the Wess-Zumino-like term in the D7 
brane probe action implies the existence of previously overlooked new phases. 

In the case 1), we have found a "chiral helix" phase in which the U(1)a angle of D7 
brane embedding increases monotonically along the direction of the magnetic field. We 
consider this as a holographic realization of the chiral spiral phase in QCD. We have 
found an axial current propagating in this phase, corresponding to the chiral separation 
effect. Previously, the existence of the chiral magnetic current at finite axial chemical 
potential has been established within the same D3/D7 model [3]. It was argued in Ref. 
[38| 139] that the coupling between the axial and vector charge oscillations induced by 
the anomaly should lead to the emergence of a gapless excitation - the chiral magnetic 
wave. It would be interesting to establish the presence of this excitation and to study its 
properties within the D3/D7 holographic model. 

In the case 2), we have identified the "spontaneous rotation" phase in which the D7 
brane spontaneously begins to rotate, so that the U(1)a angle chang function of 

time. This phase is a geometrical realization of a phase with non-zero chiral chemical 
potential - indeed, the parallel external electric and magnetic fields generate chirality 
through the triangle anomaly. In the D3/D7 model the U(1)a symmetry is shared by 
other adjoint matter fields, so the U(1)a charges in the hypermultiplet can be lost to the 
adjoint sector. This loss of chirality eventually leads to a stationary rotation frequency 
of the D7 brane, corresponding to some limiting value of the chiral chemical potential 
that can be maintained in the system. In real QCD plasma, the axial charge can be lost 
due to the sphaleron transitions; it would be interesting to establish whether or not there 
exists a corresponding limiting value of the chiral chemical potential similar to the one 
we observe in the D3/D7 model. In general, our study indicates an important role played 
by the anomaly in the phase diagram of gauge theories. 
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